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Abstract — This work considers worst-case utility maximization 
(WCUM) problem for a downlink wireless system where a multi- 
antenna base station communicates with multiple single-antenna 
users. Specifically, we jointly design transmit covariance matrices 
for each user to robustly maximize the worst-case (i.e., minimum) 
system utility function under channel estimation errors bounded 
within a spherical region. This problem has been shown to 
be NP-hard, and so any algorithms for finding the optimal 
solution may suffer from prohibitively high complexity. In view 
of this, we seek an efficient and more accurate suboptimal 
solution for the WCUM problem. A low-complexity iterative 
WCUM algorithm is proposed for this nonconvex problem by 
solving two convex problems alternatively. We also show the 
convergence of the proposed algorithm, and prove its Pareto 
optimality to the WCUM problem. Some simulation results are 
presented to demonstrate its substantial performance gain and 
higher computational efficiency over existing algorithms. 

I. Introduction 

Linear transmit precoding has been recognized as an impor- 
tant technique for capacity improvement and low implemen- 
tation complexity. Considering a transmit design for system 
utility maximization in a single-cell multiuser multiple-input 
single-output (MISO) wireless system, several works have 
focused on finding the optimal beamforming solution for the 
problem [l]-[3]. However, since the problem is NP-hard in 
general [4], the computation complexity of any algorithms 
for finding the optimal solution can be prohibitively high, 
rendering the convergence rate very slow. Therefore, those 
algorithms may be infeasible for real-time implementation. In 
light of this, considerable attention has been shifted towards 
finding a more accurate and computationally efficient solution 
for the system utility maximization problem. Assuming that 
the base station (BS) can perfectly acquire the channel state 
information (CSI) of the users, some efficient suboptimal 
algorithms have been proposed for the utility maximization 
problem [5]-[7]. The authors in [7] propose a suboptimal al- 
gorithm based on the idea of successive convex approximation 
(SCA) for a weighted sum rate maximization problem, with 
numerical results showing that their algorithm outperforms 
existing ones. 

However, in practical situations, it is inevitable to have 
channel estimation errors at the BS due to finite training 
resource (e.g., power and signal length) or limited feedback 
bandwidth [8]. If one uses inaccurate CSI directly to design 



the transmit precoders, then the system performance may be 
degraded seriously. In view of this, we design a transmit 
precoder by taking CSI errors into account. Also, in a slow 
fading channel, under a strict constraint on quality of service 
(QoS), the system must be designed for the worst-case scenario 
[9]. In this paper, the CSI errors are assumed be bounded 
spherically, and the transmit covariance matrices for all the 
users are designed to maximize the worst-case (i.e., minimum) 
system utility function against any possible CSI error, subject 
to a total transmit power constraint. However, the worst-case 
utility maximization (WCUM) problem is nonconvex and is 
generally hard to efficiently solve. Therefore, very few efficient 
algorithms have been reported for the WCUM problem. 

One approach to developing a suboptimal solution for the 
WCUM problem is based on SCA. Although the SCA-based 
algorithm presented in [7] relies on the assumption that the BS 
has perfect CSI of the users, the SCA-based concept can be 
easily extended for handling the WCUM problem. Specifically, 
we can apply a conservative convex approximation to the 
WCUM problem, and then iteratively update the associated 
parameters to improve the performance. Since a conservative 
approximation is used in SCA-based algorithms, performance 
of such algorithm would be degraded. Hence in this paper, we 
propose a low-complexity algorithm for the WCUM problem 
without any conservative approximations. In the proposed 
WCUM algorithm, two convex optimization problems are 
solved alternatively. The proposed algorithm is guaranteed to 
converge, and the limit point is Pareto optimal to the WCUM 
problem. Some simulation results are presented to show that 
the proposed WCUM algorithm significantly outperforms the 
SCA-based algorithm, and the former has lower computational 
complexity than the later. 

II. System Model and Problem Statement 

A. System Model 

We consider a single-cell multiuser multiple-input single- 
output (MISO) wireless system, where a BS equipped with 
N t antennas communicates with K single-antenna users. The 
transmitted signal at the BS is given by 



x(t) 



K 

£ 

fc=i 



Xk(t), 



(1) 



where x k (t) G C N * denotes the information-bearing signal for 
the fcth user. With (1), the received signal at the fcth user can 
be represented as 

y k (t)=h%x(t)+n k (t), (2) 

where hk G C ' denotes the channel vector between the BS 
and the fcth user, and n k (t) is the additive noise at the fcth 
user with power a\ > 0. 

Assuming that x k (t) is complex Gaussian distributed with 
zero mean and covariance matrix Qk h (positive semidefi- 
nite (PSD) matrix), i.e., x k (t) ~ CAf(0, Qk), and considering 
single-user detection, the achievable rate of the fcth user can 
be represented as (in bits/sec/Hz): 

flfc ({q«}£x, h k ) = log 2 ( 1 + A (3) 

The goal of the transmit precoding is to design the transmit 
covariance matrices {Qk} k =i sucri mat a s y stem utility func- 
tion U(Ri, . . . , Rk) can be maximized, as formulated as the 
following optimization problem: 

max C/({i? fc ({QjEi,^)}f=i) (4a) 

Qk EM * , 
k=l,...,K 

K 

s.t. ^Tr(Q fe )<P, Q 1} ...,Q K hO, (4b) 

k=l 

where Tr(-) denotes the trace of a matrix, and P > is a 
preset maximum allowed total transmit power. We assume 
the system utility function U(Ri,...,Rk) to be strictly 
increasing and concave with respect to Rk, for k = 1, . . . , K, 
as satisfied by many practical system performance measures 
[4], e.g., sum-rate utility U{R U ...,R K ) = (1/K) £f =1 R k . 
Although problem (4) is not convex in general due to the 
nonconcave utility function with respect to {Qk}k—i> ^ ^ a ^ s 
in the class of monotonic optimization problems [10], and the 
optimal transmit covariance matrices can be obtained [l]-[3], 
in spite of high computational complexity. 

B. Problem Statement 

In practice, the BS cannot perfectly acquire the CSI from 
users due to finite training power or limited feedback band- 
width [8]. In this work, the true channel is modeled as 

h k = hk + e k , (5) 

where hk G C N * denotes the channel estimate of hk at the BS, 
and efe G C * is the corresponding CSI error vector which 
is assumed to lie in a norm ball with radius rk > 0, i.e., 
\\ e k\\ < ^k, where || • || denotes the vector 2-norm. 

Our goal is to jointly design the transmit covariance matrices 
{Qk\k=i sucn mat me worst-case (i.e., minimum) system util- 
ity over the CSI errors is maximized, subject to a total transmit 
power constraint. Mathematically, the WCUM problem can be 



formulated as 

max min U {{R k {{Qi]tiM + e k )}^ =1 ) (6a) 

Q k en N t, IK||<r fe , 

k=l,...,K k=l K 

K 

s.t. ^Tr(Q fe )<P, Qi,...,Q K hO, (6b) 

fc=i 

Equivalently, problem (6) can be reformulated as 

max U(t\,...,tK) (7a) 

k=l,...,K 

Bt (h k +e k ) H Qk{h k + e k ) 

E£ fc (fc* + e k ) H Qe(hk + e k ) + a 2 k ~ *' 
V|[e fe || <r k , k=l,...,K, (7b) 

K 

Tr (Qfe) < p ' Qi,- -,QKhO, (7c) 

fe=i 

where we have made the change of variables 

R k ({Qi}? =1 ,h k + e k ) = log 2 (l + t k ) (8) 

in the utility function. Problem (6), or equivalently problem 
(7), provides a maximum lower bound on the system utility 
function over the CSI errors, given that {hk}^ =1 is known at 
the BS. However, problem (7) is hard to solve due to the in- 
finitely many nonconvex constraints in (7b). Before presenting 
the algorithm for efficiently and effectively handling problem 
(7), let us present a simulation example to demonstrate the 
essentiality of the WCUM design. 

Example: Consider sum-rate utility for problem (7) and 
naive-CSI-based design. For naive-CSI-based design, the 
transmit covariance matrices are obtained by solving the con- 
ventional perfect-CSI-based problem (4), where the channel 
estimates {hk}^ =1 are used as if they were true channels. For 
each realization of the channel estimates {hk}^ =1 generated 
with complex Gaussian distribution, the optimal transmit co- 
variance matrices for problem (7) and naive-CSI-based designs 
are obtained by applying branch-reduce-and-bound (BRB) al- 
gorithm reported in [3]. With the obtained transmit covariance 
matrices and the presumed channels {hk} 1 ^^, the worst-case 
sum rates were determined to be the minimum value of Rk 
computed by (3) over 10 6 true channels {hk = hk + e.k}^ =l , 
where the simulated CSI errors {ek}^ =1 are randomly and in- 
dependently generated satisfying ||ei.|| < r for fc = 1, . . . , K. 
Figure 1 displays some simulation results of the achieved 
worst-case sum rate versus r for N t = K = 2, P = 10 dB, 
and erf = • • • = a 2 K = 0.01, where each result was obtained 
by averaging over 100 realizations of the channel estimates 
{hk}k=i- F rom this figure, one can see that improper transmit 
covariance matrices can yield a very low sum rate in the worst 
case, especially for large CSI error radius. ■ 

Since problem (7) is NP-hard in general [4], the compu- 
tational complexity of any algorithms for finding the optimal 
solution of problem (7) can be prohibitively high, and thus 
it is infeasible for real-time implementation [l]-[3]. Next, 



WCUM problem (7) 
Naive-CSI-based design 




Fig. 1: Achieved worst-case sum rate vs. CSI error radius r. 

we concentrate on the proposed low-complexity algorithm for 
finding a more accurate suboptimal solution for problem (7). 

III. A Low-Complexity WCUM Algorithm 

In this section, we first present the proposed iterative 
algorithm for problem (7), and then prove its convergence. 

A. Proposed Algorithm 

To proceed, let us apply S -lemma [11] to constraint (7b), 
and then problem (7) can be equivalently reformulated as (see, 
e.g., [3] for details): 



max U(ti, 

k=l,....K 



,t K ) 



s.t. *fc(tfc,A fc) {Qi} 



K \ 
i=i 



y o. k = i, 



,K, 



K 

E 

fc=i 



Tr(Q fe ) < P, 



A fc >0, Q k h0, k = l,...,K, 



(9a) 

(9b) 
(9c) 
(9d) 



where Ai, . . . , Xk € I 



AfcJjV t 





are the introduced slack variables, and 

H 



J-Nt 



-X k r\ - t k a\ 



1,...,K. (10) 



Although problem (9) is still not convex due to nonconvex 
constraint (9b), the problem has a more tractable form than 
problem (7). To develop the proposed algorithm, we need the 
following lemma. 

Lemma 1 Each eigenvalue of the matrix (t k , X k) {Qi}|£i) 
defined in (10) decreases with t k , for k = 1, . . . , K. 

Proof: By letting S kjj G R, j = 1, . . . , N t + 1, be the 
eigenvalues of the matrix & k (t k , X k , {Qi}|Li) and denoting 



have 



■v/V t + l 



as the associated unit-norm eigenvector, we 







v kij Y kik v ktj + X,v!', 



v k , 



\ e^k 

where [a]j denotes jth entry of a vector a, and 

H 



< k ,e 



^0, £ = !,. 



(11) 



(12) 



Since Yfe^ is a PSD matrix, from (11), one can easily show that 
the eigenvalue 8 k ,j decreases with t k , for j = 1, . . . , N t + 1. 
This proof is thus complete. ■ 
Problem (9) aims to maximize t k (because the utility func- 
tion is strictly increasing in t k ). According to Lemma 1, one 
can maximize t k by maximizing the minimum eigenvalue of 
the matrix $ k (t k , X k , {Qi}fLi) i n (9b) (which only involves 
t k ). Also, it can be inferred that the minimum eignevalue of 
&k(t kl X k , {Qi}fLi) must be zero as the optimal solution of 
problem (9) is achieved. Based on these facts, let us consider 
the following iterative approach for dealing with problem (9). 
At the mm iteration, ti,...,tic are updated by solving 



max 

fc=l K 



U(tx,...,t K ) 



s.t. * fe (44 m \{Qr}£i)^0, Vfc 



(13a) 



(13b) 



Let us denote the obtained optimal utility value in problem 



,K, are obtained by solving 



(13) as U^dQr^Xr'}^), in which Q k 



(m) 



and A 



(m) 



max 

k=l,...,K 



fc=i 



K 

E 

/c=i 



Tr(Q fc ) < P, 
A fc > Q fe >r 0, z fc > 0, fc 



(14a) 

»k-Wi>:0, Vfc, (14b) 
(14c) 

1,.,.,K. (14d) 



The obtained optimal value of in problem (14) is denoted 



as * (m) ({4 m_1) }f =1 ), in which t^ 1 ', k = 1,...,K, 
represent the solution obtained by solving problem (13) from 
the (to — l)th iteration. Problems (13) and (14) are convex, 
and thus can be efficiently solved. Specifically, since the 
design variables t\, . . . ,tx in the objective function and in the 
constraints of problem (13) are decoupled, and the objective 
function is strictly increasing in t k , for fc = 1,...,K, by 
Lemma 1, the optimal {t k }^ =1 can be separately obtained by 
simple bisection search; i.e., find a t k such that the minimum 
eigenvalue of $> k (t k , Ajj, , {Q,- }|Li) in (13b) is equal to 
zero. The obtained WCUM algorithm for problem (7), or 
equivalently problem (9), is summarized in Algorithm 1. 



Algorithm 1 Proposed WCUM algorithm for problem (7). 

1: Obtain a feasible point {ij^}|Li according to (17); set 
solution accuracy e > 0; and set iteration index m = 0. 

2: repeat 

3: Update m := m + 1. 

4: Obtain {Q { ™\ X k m) } k=1 by solving problem (14). 

5: Obtain {t { k m) }^ =l by solving problem (13). 

6: until the predefined stopping criterion is met, e.g., | £A m ) 

£f(m-l)| < £ _ 



An initial feasible point {t k } k= i f° r solving problem (14) 
in the first iteration can be obtained by finding a feasible point 
of problem (7). Let the transmit covariance matrices be of 
rank one, i.e., Q k = w k w k for k = 1, . . . , K, where w k G 
C k = 1, . . . , K, can be any arbitrary vectors such that the 
power constraint in (7c) is satisfied, i.e., EfcLi ll^fcll 2 — P- 
Therefore, from (7b), a feasible point {t k } k=1 is given by 



f (0)_ 

l k — 



minj| efc ||< rfc \(h k + e k ) H w k \ 2 
max IK||<r fc E^fe \(hk + e k ) H w t \ 2 + <r\ 

([\h^ Wk \~r k \\w k \\]+f 
Ef^(l^l+^ll^ll) 2 +^, 



k = 1, 



where [a] + = max{a, 0}. Letting 




w k = 

in (15) gives rise to 

f( o) _ j([\\hf\\-r k ] 



k = l,...,K, 



,K, (15) 



(16) 



t. 



p_ ( |hf fe< 
L-W K \ \\h t \\ 



'■/,• 



1,...,K. (17) 



B. Convergence of Algorithm 1 

In the following, we will show the convergence of Algo- 
rithm 1 and the limit point to be Pareto optimal to problem (7). 
To this end, we need the following lemma: 

Lemma 2 If the objective value of problem (14) obtained in 
the mth iteration is positive, i.e., <3>( m ) > 0, then the system 
utility value can be increased by solving problem (13), i.e., 
Jj( m ) > LA™ -1 ), and the limit values z* , . . . , z* K obtained by 
Algorithm 1 must be zero. 



Proof: Since \]/( m ) = J2 k =i z k I s positive, let us assume 
> for some k without loss of generality. Therefore, 



Am) 



from (14b), we have 

According to Lemma 1, we can always find a value, say t 



(18) 

(m) 



such that t 



(m) 



> t 



(m-1) 
fc 



is feasible to problem (13). As a 



result, we have UW > L^™" 1 ) since the objective function 
of problem (13) is strictly increasing in t k . 



Next, let us show that as Algorithm 1 converges, the limit 
values Zi,...,Zj( are all zero by contradiction. If z£ ^ 0, 
then we can further increase the objective function of problem 
(13), which contradicts the premise that the algorithm has 
converged. This proof is thus complete. ■ 

Since LA m ) in Algorithm 1 is monotonically increasing 
in the iteration number m [by Lemma 2] and its value is 
bounded above due to finite total transmit power P, we can 
conclude that Algorithm 1 must converge. In the following 
proposition, the Pareto optimality of Algorithm 1 to problem 
(7) is established. 

Proposition 1 The limit point {Q k , ^ k } k= i obtained by 
Algorithm 1 is Pareto optimal to problem (7). 

Proof: To show that the limit point {Q k ,t k , X k } k=1 is 
Pareto optimal to problem (7), or equivalently problem (9), we 
first show that the limit point is a feasible point of problem 
(9). Since the point {Q k , X k , z^, tt}^ = i is feasible to problems 
(13) and (14), we have 



* fe M,{Q*}f=i ) >=o, k = i 



K 

E 

fc=i 



Tr(Q*) < P, 



(19a) 
(19b) 
(19c) 



K > 0, Ql hO, k = l, 

where we use the fact z\ — ■ ■ ■ — z* K = in (19) according 
to Lemma 2. Comparing (19) with the feasible set of problem 
(9), one can conclude that the limit point {Q k , t k } k =i 
generated by Algorithm 1 is feasible to problem (9). 

Now, let us show that the limit point {Q k ,t%, \ k } k=1 
is Pareto optimal to problem (9) by contradiction. Suppose 
that the limit point is not Pareto optimal to problem (9). 
Then, according to the definition of Pareto optimality in 
[10], there exists a feasible solution to problem (9) such that 
U(h, ...,t K ) > U{t\, t* K ), where i k > t% and l t > t\ 
for £ 7^ fc. That is to say, the following constraint set is feasible 
[by Lemma 1]: 



* t (t* e ,\ t ,{Qi}? =1 ) h 0, W^fc, 

K 

^Tr(Q fc ) <P, 

k=l 

A fe > 0, Q k hO, k = 1, . 



, K. 



(20a) 
(20b) 

(20c) 

(20d) 



The feasibility of constraint (20) implies that there exists a 
solution {\,Qi}fLi sucn that the minimum eigenvalue of 
the matrix &k{t%, Afe, {Qi}|Li) is positive, rendering z\ > 
in problem (14), which contradicts with the fact of z£ = 
by Lemma 2. Therefore, we have completed the proof that 
the limit point yielded by Algorithm 1 is Pareto optimal to 
problem (7). ■ 

IV. Simulation Results 

We consider the wireless system as described in Section II 
with Nt = 4 transmit antennas at the BS and K = 2 single- 
antenna users, the total transmit power P = 10 dB, and the 



users' noise powers <j\ = erf = 0.01. For simplicity, the CSI 
error radii of all users are assumed to be identical, i.e., r\ = 
■ ■ ■ = Tk = t. In each simulation trial, the presumed channels 
{hk}ff = x are randomly and independently generated according 
to the standard complex Gaussian distribution. 

Considering sum-rate utility for problem (7), we compare 
the worst-case sum rate performances of the proposed WCUM 
algorithm (Algorithm 1), SCA-based algorithm [7] and the 
optimal (i.e., maximum) worst-case sum rate obtained by 
BRB algorithm (which is a brute force approach) reported 
in [3]. The solution accuracy for Algorithm 1 and SCA-based 
algorithm is set to 10~ 3 , i.e., e = 10~ 3 , and the gap tolerance 
between the upper and lower bounds for the BRB algorithm 
is set to 0.1 as in [3]. The involved convex problems in the 
algorithms under test are solved using CVX [12]. Figure 2 
shows the average sum rate versus CSI error radius r, where 
the sum rates obtained from the three algorithms are averaged 
over 100 realizations of the presumed channels {hk}^ =1 - 
From this figure, one can see that the proposed algorithm 
performs much better than SCA-based algorithm, with the 
performance gap about 1.8 bits/sec/Hz, but worse than BRB 
algorithm, with the performance gap between 0.8 bits/sec/Hz 
for r — 0.1 and 0.2 bits/sec/Hz for r = 0.2. Note that the 
sum rate yielded by the proposed algorithm is closer to the 
optimal (i.e., maximum) sum rate as the error bound r grows, 
and that the performance gap between SCA-based algorithm 
and BRB algorithm is as high as 2.6 bits/sec/Hz for r = 0.1 
and 2 bits/sec/Hz for r = 0.2. 

By our simulation experiences, we found that the proposed 
algorithm converges much faster than BRB algorithm and 
is more computationally efficient than SCA-based algorithm. 
As an illustration, the average computation times of BRB 
algorithm, SCA-based algorithm, and the proposed algorithm 
for obtaining the results shown in Fig. 2 for r = 0.1 using a 
desktop PC with 3GHz CPU and 8GB RAM, are 550.9 sees, 
123.9 sees, and 29.6 sees, respectively. The fast convergence 
rate of the proposed algorithm is due to the use of simple 
bisection method for solving problem (13), and due to the 
smaller size of problem (14) compared with those problems 
involved in BRB and SCA-based algorithms. Let us emphasize 
that for large values of N t and K, BRB algorithm will not be 
applicable due to extraordinarily high complexity. However, 
the detailed computation complexity analysis of Algorithm 1 
is omitted here due to the space limitation. 

V. Conclusion 

We have presented a low computation complexity algorithm 
[see Algorithm 1] for finding a more accurate suboptimal so- 
lution for the WCUM problem in (6). The proposed algorithm 
has been proved to converge, and the limit point is Pareto 
optimal to the problem [see Proposition 1]. The presented 
simulation results have demonstrated that the proposed algo- 
rithm performs much better than SCA-based algorithm and 
has higher computational efficiency over both BRB and SCA- 
based algorithms. 
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Fig. 2: Performance comparison of the proposed WCUM 
algorithm and two existing algorithms in terms of average sum 
rate versus CSI error radius r. 
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